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Abstract: We show that the Bianchi I Einstein field equations in a perfect fluid scalar 
field cosmology are equivalent to a linear Schrodinger equation. This is achieved through 
a special case of the recent FLRW Schrodinger-type formulation, and provides an alternate 
method of obtaining exact solutions of the Bianchi I equations. 
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1. Introduction 

Recently, a correspondence was established between solutions of Einstein's field equations 
in a Friedmann-Lemaitre-Robertson- Walker (FLRW) universe and solutions of a particular 
nonlinear Schrodinger-type differential equation. That is, given a solution of the latter, a 
solution of the former can be constructed via the prescription given in [3] (and vice versa). 
Further motivated by a connection between Bianchi I and FLRW cosmologies seen in a 
paper by James E. Lidsey [6], an analogous linear Schrodinger formulation is demonstrated 
here for the anisotropic Bianchi I universe. The author would like to extend many thanks 
to Floyd L. Williams for his valued advice on the results presented here. 

2. Einstein Equations 

Consider the Einstein field equations = K^Gij for a perfect fiuid Bianchi I universe with 
scalar field 0, potential V and metric ds'^ = -dt^ + xltfdx'^ + Y{t)'^dy'^ + Z{tfdz^. We 
will consider the case where the energy density and pressure are given solely by a scalar field 
contribution, i.e. no matter contribution. That is, p — 0^/2 + y o and p ~ (/'^/2 ~ V o cj). 
For a vanishing cosmological constant the equations take the form 
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where K'^ = SttG and G is Newton's constant. 

The fluid conservation equation can be derived from these equations and is 

p + e{p+p) = 

where 

-(144) 

is the expansion/contraction of vohime. By definitions of p and p, this equation reduces to 
the Klein- Gordon equation of motion 

^ + 0<j) + V' o(j) = 0. 

Note that for 7 = 2(/>^/[0^ + 2{V o (f))] one has the equation of state p = (7 — 

3. Description of The Correspondence {X, Y, Z, (p, V) < — > u 

Similar to the formulation in Lidsey [6], we first describe how a solution to the Bianchi 
I equations (i)-(iv) can be used to derive a solution to the FLRW equations. Since [3] 
provides the FLRW-Schrodinger connection, this will motivate the Schrodinger-Bianchi I 
correspondence. 

We begin by defining the quantities 

X Y X Z Y Z 

Vi^^-y, m^j^-^, Vs^y-^. (3.1) 

Computing ■^[{i) — {ii) — {in) — {iv)] and using the definitions above, one can verify Ray- 
chaudhuri's equation 

e + 2ii''-^e^ + ^{2p + p) = Q (3.2) 

where p is the shear scalar given by /x^ = | {t]\ +r]\ + rj^) + |^^. Also, rewriting equation 
(i) using the above notation, 

le'-p'' = K^p. (3.3) 

Remarkably, equations (3.2) and (3.3) are a special case of the FLRW field equations 
(see [3]) with the following substitutions: 

n = 6, fc = 0, a{t) = {XYZf\ D = ^^^^ {^f + rj^ + 4) . (3.4) 

Note that [3] requires D to be constant. By (i)-(iv) in (2.1) and the following lemma, one 
can easily show that each of the products XYZrji for i = 1, 2, 3 is a constant function of t. 

Lemma. For arbitrary functions X{t), Y{t), Z{t) > and f{t), 

f + 6f = fXYZ is a constant function 

for e as in (2.2) in terms ofX, Y, Z. 
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Therefore, given a quintet {X,Y, Z,(j),V) solution to (i)-(iv) in (2.1), one can construct a 

solution to the FLRW field equations by (3.4). By the converse of Theorem 1 from [3], one 
can then construct a solution to the time-independent linear Schrodinger equation 



u"{x) + [E- P{x)]u{x) = 0. (3.5) 



with constant energy E and potential P{x). We will see that the Schrodinger solutions 
derived in this way will always be such that _B < 0. 

In this paper, we show the direct correspondence (X, Y, Z, (j>, V) < — > u between solutions 
{X,Y,Z,(j),V) of (i)-(iv) in (2.1) and solutions u of (3.5). This correspondence provides an 
alternate method of solving Bianchi I field equations. 

With the above notation in place, we can now state the main theorem: 

Theorem. Let u{x) be a solution of equation (3.5), given E < and P{x). Then 
a solution {X,Y, Z, (f),V) of the Einstein equations (i)-(iv) in (2.1) can be constructed as 
follows. First choose functions <j(t), tp{x) such that 

a(t) = u{a{t)), ^'(xf = ^P(x) (3.6) 
and also constants ci, C2 such that 

cl + ciC2 + cl = — —. (3.7) 

JVext define the functions 

R(t) = u(a(t))-^/^ (3.8) 

and 

"(i) = |^(i), P{t) = ^a(t), ^(t) = -a{t)-l3it)- (3-9) 
Then the following quintet solves Einstein's field equations (i)-(iv): 

X(i) = i?(t)e«W, F(i) =i?(i)e^W, Z(t) = i?(i)eTW, (3.10) 

</.(*) =V(a(i)), V=^[(u'r+u'[E-P]\oi,-\ (3.11) 
Here, in fact, (X, Y, Z, (j), V) will also satisfy the equations 
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Vim = ^ {(>'+(>)■ (3-13) 

Conversely, let (X,Y,Z,(j),V) be a solution of equations (i)-(iv) in (2.1), with p and p 
as before. Similar to (3.6), choose some solution a{t) of the equation 



m = (3.14) 



Then equation (3.5) is satisfied for 



E = -^X^Y^Z\ii+ii+4), (3.15) 
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Note that in (3.15), E < and is constant by the same argument stated for D above 
the lemma. The theorem therefore provides a concrete correspondence (X, Y, Z, cf), V) ^ u 
between solutions {X,Y, Z,(j),V) of the field equations (i)-(iv) and solutions u of the linear 
Schrodinger equation (3.5). 
Remarks. 

1. The case P{x) = 0: 



In most examples P{x) is nonzero. However, if P{x) = then the theorem must be 
stated carefully, as ipix) is a constant function by (3.6) and has no inverse. Therefore 
the expression for V in (3.11) has no moaning. In this case, we will show that the 
right-hand side of (3.13) is a constant function that will serve as our new definition for 
V. By (3.8)-(3.10), uo(j= 1/{XYZ). Differentiating and using (3.6), {u'o(7)(uoa) = 
—9/{XYZ). That is, u' o a = —9. Differentiating again, {u" o a)(u o a) = —6. Using 
these in (3.5), composed with a and multiplied by u o cr, 

^ E , ^ 

= 0. (3.16) 



^2^2 ^2 



Therefore (3.12) is still valid, of course, with the left side equal to zero since <f) is 
constant in this case by (3.11). Differentiating (3.16), 

=0. (3.17) 



Now, to show that (3.13) is constant, differentiate its right side to get 

d 



= by (3.16). 



Therefore in the case P{x) = when the definition of V in (3.11) no longer has 
meaning, we define V{x) = Vq = {6^ + 6)/{3K^) and (l){t) = any constant; and we note 
that equations (3.12) and (3.13) still hold in this special case. 

2. Equations (3.12) and (3.13) imply (i)-(iv) in (2.1) under a condition: 



By the comment preceding the lemma above, any solution to the equations (i)-(iv) will 
have the property that XYZr]i is constant for rjj as in (3.1) and i = 1,2,3. Suppose 

we are given a priori (positive) functions X, Y, Z with this property. Differentiating 
each of XYZrji and setting equal to zero shows exactly that the left sides of (ii)-(iv) 
are equal to each other. Next, any three positive functions can be reparametrized as 

in (3.10) for R{t) = {XYZf'^, a = i In , ^ = 5 In and 7 = i In (fl) . 

Using these formulas and the constant quantities XY Zrji^ one can easily compute that 
each of d,/3 and 7 are scalar multiples of XjXYZ^ therefore establishing (3.9) for a as 
in (3.14). Finally, under the condition that Ci,C2 satisfy (3.7) for the constant < 
given by X^Y,Z as in (3.15), one can show that equations (3.12) and (3.13) indeed 
compute and V solving (i)-(iv) in terms of the given X, Y, Z. 
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4. Examples 



As an illustration of the theorem, take the solution u{x) = Ae~^^~^^ — i?e^^~^^ for A,B>0 
to the equation (3.5) with < and P{x) = 0. Solving the differential equation (3.6) for a 
using Mathematica, we obtain 



a{t) 



:ln 



\WL\^\^J-ABE{t- Cq) 



(4.18) 



for integration constant cq. Also by (3.6), V(a;) = V'o = any constant. Then by (3.8) 

1/3 



-1/3 



sinh[2V^ARE(t - co)]^ 



(4.19) 



Further let ci,C2 be any constants such that (3.7) holds given the constant choice E. We 
form X, Y, Z, <f) according to (3.9)-(3.11) and obtain 



X = (^^-^sinh[2\/=ABE(i-co) 
Y = (^ ^Jj^ smh[2V- ABE{t - cq) 
Z = (^ ^Jj^ smh[2\/ -ABE{t - cq) 
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ci/(2v^) 



— tanh[y/-ABE{t - cq)] 
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— tanh[V-ABE{t - cq)] 



tanh.[\/ -ABE{t - cq)] 



C2/(2v^) 



-(ci+C2)/(2v^) 



for t > Co and ^ = Vo- Since P = 0, does not exist (see Remark 1) and we use (3.13) 
for the definition of V and obtain the constant V = Vo = —AABE/ {2>K^). The reader may 
compare this solution with a similar one in [1]. 

As another example, we will begin with the same assumptions on E, P and will obtain 
quite a different Einstein solution. That is, again let S < and P{x) = 0, but take solution 
u{x) = Aer^/^^ to (3.5) with A > 0. Solving the differential equations in (3.6), wo obtain 
a(£) = In [A^/—E{t — co)] l\/~-E for t > co and i!{x^ = -00 = any constant (therefore we 
will also have = ^/jq by (3.11)). By (3.8), R(t) = {V^{t - cq))^/^. Letting ci,C2 be any 
solution to (3.7), finally we compute X,Y, Z to be 
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Y 

Z 



= A^^/(2v^) (x/=:B(t-co)) 

^ ^-(cx+c.)/(2y=ij) (y^(^t _ Co)) 



ci/(2v^)+l/3 



(4.21) 



C2/(2\/=E)-|-l/3 



-(ci+C2)/(2v^)+l/3 



by (3.9)-(3.10) for t > cq. Again, since P = 0, ip~^ does not exist and we use (3.13) as the 
definition for V and obtain ^ = 0. That is, this solution is vacuum. 

To further demonstrate the utility of the theorem, we will take a trivial non-physical 
solution of (i)-(iv), and first use the converse theorem to map to a solution of (3.5). We will 
then apply the theorem a second time and map back to a solution of (i)-(iv) and will have 
produced a physically acceptable example. Begin by considering the vacuum = V = 0) 
solution 

X = Roe°"', Y = Roe'^°, Z = Roe^° (4.22) 
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for constants Rq > 0,ao, /?o, 7o such that ao + /3o + 7o = 0. By (3.1) and (3.15), E = P = 
and u{x) mq = (I/Rq). Clearly this is a solution to the linear Schrodinger equation (3.5). 
Note that since u is constant and (p is zero, we did not need to compute a. Now to map back 
to a solution of (i)-(iv), we solve (3.6) and use (3.11) so that a{t) = {1/R^)t and ip = (p = 
any constant. Now by (3.8)-(3.10), 



(4.23) 



Again by (3.13), V = 0. 

As a final example, wc take u{x) = {l/x)e^^ E < and P{x) = (2/a;^) + E'^x'^ for 
a; > 0. This Schrodinger solution was found using the techniques in [5]. Solving (3.6), we 

obtain a{t) - 



ln[—E{t — Co)] for integration constant cq and 



ipix) 



1 



y/6K 



\/2 + E'^x^ + V2hi 



2 + V4 + 2E^3 



(4.24) 



Graphing this function for a few values of E indicates that the inverse exists, and we will 
denote it by Calculating X, Y, Z, cp, V according to the theorem. 



X = 

Y = 
Z = 

V = 



(^l^t-co)^/-2E\n[-Eit-co)]) e'='/V'-[-^(t--o)]/(-2E) 
((i - Co)V-2Eln[-E{t-Co)]) ^'^ e^^/ViW-iJ(t-co)]/(-2iJ) 
((i - co)^J-2E\n[-E{t-co)]) e-(^i+'==)/Vi"[-^(*--o)]/(-2E) 



(4.25) 
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(1 + Ex^) 



^l + 21n2[-i;(i-Co)]+ln 



■ ln[-£;(t - Co)] 



for ci, C2 satisfying (3.7) and t > cq. 



E + EJI + 2 In^ [-£;(* - Co) 
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